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We apply two methods, i.e., the Gaussian processes and the non-parametric smoothing procedure, 
to reconstruct the Hubble parameter H{z) as a function of redshift from 15 measurements of the 
expansion rate obtained from age estimates of passively evolving galaxies. These reconstructions 
enable us to derive the luminosity distance to a certain redshift z, calibrate the light-curve fitting 
parameters accounting for the (unknown) intrinsic magnitude of type la supernova (SNe la) and 
construct cosmological model-independent Hubble diagrams of SNe la. In order to test the com¬ 
patibility between the reconstructed functions of H{z), we perform a statistical analysis considering 
the latest SNe la sample, the so-called JLA compilation. We find that, for the Gaussian processes, 
the reconstructed functions of Hubble parameter versus redshift, and thus the following analysis on 
SNe la calibrations and cosmological implications, are sensitive to prior mean functions. However, 
for the non-parametric smoothing method, the reconstructed functions are not dependent on initial 
guess models, and consistently require high values of Ho, which are in excellent agreement with 
recent measurements of this quantity from Gepheids and other local distance indicators. 

PACS numbers: 


I. INTRODUCTION 


The inherent relation between the peak luminosity of type la supernovae (SNe la) and the speed of 
luminosity evolution after maximum light (known as the Phillips relationship) [l|, makes possible to 
standardize these events as a distance indicator and measure the geometry and dynamics of the universe. 
As well known, several years after the discovery of the Phillips relationship, observations of some dozens 
of distant SNe la led to the discovery of the cosmic acceleration M- In Einstein’s general relativiW, 
such behavior implies either the existence of a new field, the so-called dark energy (see Refs. 5, 0 
for recent reviews), or that the matter content of the universe is subject to dissipative processes 3-[3- 
The mysterious cause of the current cosmic acceleration can also be attributed to a modification of the 
standard theory of gravity itself on cosmologically relevant physical scales [iMl. 

In the past decade, several groups have put a lot of effort into enlarging the sample size of well-measured 
SNe la events. At the same time, improvements in precision and assessments of systematic errors have 
also been accomplished. Recently, several SNe la samples containing a large number of SNe la events 
with high quality have been released, such as the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II 
SN Surve y) the Union2 and Union2.I SN la released by the Supernova Cosmological Project 

(SCP) the first three years of Supernova Legacy Survey (SNLS3) [l3)[3> the joint light-curve 

analysis (JLA) of SDSS-II and SNLS3 [ll, [2^. 

The distance estimation from SNe la data is based on the empirical observation that these events form 
a homogeneous class whose remaining variability is reasonably well captured by two parameters. One of 
them describes the time stretching of the light-curve (xi) whereas the other describes the SNe la color 
at maximum brightness (c). In the latest JLA SNe la sample [2^, which uses the SALT2 model to 
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reconstruct light-curve parameters (xi, c, and the observed peak magnitude in rest frame B band mg), 
the distance estimator assumes that SNe la with identical color, shape, and galactic environment have 
on average the same intrinsic luminosity at ah redshifts. This assumption can be quantihed by a linear 
expression, yielding a standardized distance modulus, which relates to the luminosity distance du via 


M = 51og [^]+ 25 : 


/3, Mb) = TOb ~ -^b + a x xi — P x c, 


( 1 ) 


where a and /3 are nuisance parameters which characterize the stretch-luminosity and color-luminosity 
relationships, reflecting the well-known broader-brighter and bluer-brighter relationships, respectively. 
The value of Mb is another nuisance parameter which represents the absolute magnitude of a hducial SNe 
and was found to depend on the properties of host galaxies, e.g., the host stellar masslMsteiiar), although 
the reason is not completely clear (^. [^. Here, we follow the procedure in Ref [l^ to approximately 
correct for this effect by a simple step function: 

_ J-^B: if Mstellar < IO^'^Mq. 

^ I Mg -I- Am, otherwise. 


In general, the light-curve htting parameters, a and /3, are left as free parameters being determined in 
the global ht to the Hubble diagram. This treatment results in the dependence of distance estimation 
on the cosmological model used in the analysis. Thus, implications derived from SNe la observations 
with the light-curve htting parameters determined in the global ht to the Hubble diagram are somewhat 
cosmological-model-dependent. 

Our goal in this paper is to construct a completely cosmological model-independent Hubble diagram 
of SNe la using observational data of the so-called cosmic chronometers [1^ , where the cosmic expansion 
rates M(z) are measured from age estimates of red galaxies without any assumption of cosmology. 


II. METHODOLOGY 

The expansion rate, B = a/a where a = 1/(1 -|- z), at redshifts z ^ 0 can be obtained from the 
derivative of redshift with respect to cosmic time, i.e., H(z) ~ difficult task here is 

to estimate the change in the age of the Universe as a function of redshift At. This can be done by 
calculating the age difference between two luminous red galaxies at different redshifts, as proposed in 
Ref. [15. This method is usually referred to as differential age and the passively evolving galaxies from 
which At is estimated are called cosmic chronometers. Currently, 21 measurements of H{z) based on this 
method (in the redshift range 0.070 < x < 1.965) have been obtained H^-flll- Although cosmological 
model-independent, these estimates rely on stellar population synthesis models whose influence on At, 
according to Ref. [15, becomes important at z > 1.2. In our analysis, we follow Ref. [15 and consider 
only 15 H(z) measurements up to z < 1.2 which, in practice, given the redshift distribution of the H{z) 
data, means z < 1.037. We also increase slightly (20%) the error bar of the highest-z point to account for 
the uncertainties of the stellar population synthesis models. This ensures that the evolution of the Hubble 
parameter as a function of redshift reconstructed in this paper is neither dependent on the cosmology 
nor on the stellar population model. 


A. Distances from H(z) measurements 


Recently, H{z) measurements were used to estimate distances by solving numerically the comoving 
distance integral for non-uniformly spaced data with a simple trapezoidal rule [Slj , 


dc = c 



dz' 

H{B) 


N 
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( 3 ) 
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The uncertainty associated to the 


bin is standardly propagated from the errors of H{z) data, 


s^ 



'Hi. 


"i+i 



(4) 


and then the error of the integral ([3|) from z = 0 to z = is tr^ = Y^=i = ^1- Naturally, the precision 
of this simple trapezoidal rule is sensitive to the uniformity of th e sp aced data as well as to the number 
of data points in a certain spaced range. As indicated in Ref. |32j |. the relative errors of this method 
decrease remarkably when the number of intervals averagely spaced in z = 0 — 1 increases. Here, we use 
two methods to reconstruct the evolution of the expansion rate with redshift from cosmic chronometer 
H{z) measurements, namely, Gaussian Processes (GP) and a non-parametric smoothing (NPS). This 
procedure enable us to achieve mo del-independent distance estimates by integrating the inverse of the 
reconstructed function using the approach of Ref. with a very small and uniform step of Az = z^+i — Zi 
(see Sec. 3). 


1. Gaussian processes 

As a powerful non-linear interpolating tool, the Gaussian processes allow us to reconstruct a function 
from data without assuming a model or parameterization for it. This method has been first pioneered for 
cosmology [13, [s^, and then has been widely used for several purposes, for instance, reconstructions of the 
equation of state of dark energy c osmography [Mill , null tests of the concordance model [13, [13| > 
Hq from cosmic chronometer data [40l|. and reconstructions of the distance-duality relation [4lj |. 

The reconstruction is based on a mean function with Gaussian error bands, where the function value 
at z is not independent of the function value at some other point z (especially when z and z are close to 
each other) and they are related through a covariance function fc(z, z). This covariance function depends 
on a set of hyperparameters and there is a wide range of possible candidates for it. As the function 
of Hubble parameter versus redshift is expected to be infinitely differentiable, we consider the squared 
exponential covariance function: 


fc(z,z) = cr^exp{ -(5) 

where the two hyperparameters cr/ and I are, respectively, related to typical changes in the function 
value and to the length scale one needs to move in input space to get significant change in the function 
value. In order to obtain the value of the function, the hyperparameters should be trained by maximizing 
the marginal likelihood which only depends on the locations of the observations. As the key steps 
summarized in Ref. [^ . another issue in constructing the GP is to decide a prior mean function. In 
order to achieve an unbiased reconstruction, we should choose a flat prior by taking a constant mean 
function into consideration. Moreover, it is reasonable and safe to decide a constant function when we 
do not have any prior information about the reconstructed result. As it was initiated and widely used 
in the literature, the best choice is a constant mean function because any prior model may introduce 
significant bias in the results. However, it is also necessary to check the dependence of the reconstruction 
and following implications on the selected prior mean function. Here, we do it by using different best-fit 
models, i.e., Einstein-de Sitter (E-D) model, AGDM, and wCDM as the prior mean function. Results are 
summarized in Table HI 

In this work, we reconstruct the Hubble parameter as a function of the redshift from 15 H{z) measure¬ 
ments of cosmic chronometers by using the GaPP (Gaussian Processes in Python)^ [^ . 


^ http://www.acgc.uct.ac.za/~seikel/GAPP/index.html 
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2. Non-parametric smoothing method 


In order to verify the influence of the reconstructing methods on the results, we also use the non- 
parametric procedure of Ref. (42l - l^ to reconstruct the H{z) function from cosmic chronometer data. 
The smoothing function taking into account the data errors is given by 

H%z, A) = HS{z) + N{z) X /C(z, z,), (6) 


where H^{z,A), H^{zi), H{zi), anizi) and A correspond, respectively, to the smoothed data, the initial 
guess model, the observed data, the error associated with the H{zi) data and the smoothing scale. The 
function N{z) is the normalization factor given by: 


A(z)-i 


E 


IC{z,Zi) 


(7) 


Given the arbitrariness in the choice of the kernel 1C{z,Zi), we tested a Gaussian function and the log¬ 
normal kernel of Ref. [1^ . We found a small difference between the reconstructed functions from both 
approaches. In what follows, we will adopt a Gaussian kernel. 

We perform a boot-strapping method applying iterative smoothing functions. In the first interaction, 
we subtract an initial guess model to the data in order to smooth the noise. Then we add back the initial 
guess model. In the next smoothing, we replace Hs{zi) by the previous smoothed H^{z) and calculate 
the value in each iteration. We stop the process when the variation of the value between two 
successive steps is ^ 0.002% (~ 1000*^ iteration). 

The reconstruction depends on the value of A. For example, for values of A < 0.6 the reconstructed 
function has many bumps, for A between [0.6,0.9] bumps disappear but the reconstruction at high redshift 
is strongly dependent on the A value whereas for A > 1.0 the reconstruction in the whole redshift interval 
considered is weakly dependent on the A value. However, very high values are meaningless because data 
points at high redshift should not be highly correlated with low redshift data points. In order to select 
the smoothing scale we minimize over A the cross-validation function, which is dehned by: 


i 

where H^^{xi\A) denotes the reconstructed Hubble parameter at z = Zi without taking into account the 
ponte {zi, H{zi)). For our 15 H{z) data discussed in the previous section the A value that asymptotically 
minimizes Eq. m is 1.4. 

As reported in Ref. [l^ , the initial guess model is not relevant since the reconstruction does not depend 
on it. Here, we also have tested such dependence by considering different best-fit model, namely, the E-D, 
AGDM, and wGDM as the initial guess one. Results are summarized in Table HU An important difference 
between the analysis reported here and the one performed in Ref. [1^ is that in the latter the iteration 
number of the boot-strapping process n is the free parameter relative to which a x^ function is minimized 
to define the Itr region (x^ = Xmin + !)• our analysis, in order to calculate the Ict confidence level we 
extrapolate the technique developed in for linear smoothing and consider the expression: 

where crH^(z) is the la error of the reconstruction, Vi is the smoothing factor (vi = N{z)K{z, Zt) / 
and is the estimate of the error variance given by 



.2 i:^{H{zi)-H^{zi)f 


( 10 ) 
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FIG. 1: a) Reconstructions of the Hubble expansion from cosmic chronometer data using Gaussian processes 
(blue region) with the prior mean function being a constant and non-parametric smoothing method (gray region) 
with the initial guess model being the best-fit ACDM. b) The quantity H{z)/{1 -|- a) as a function of z for both 
reconstructions. 


In order to test the validity of this extrapolation, we simulate H{z) data using the ACDM scenario as 
fiducial model with different values of cr (cTsim) and compare them with those ((Tree) calculated from Eq. 
([9]). We find that dree is — 30% smaller than dree- We, therefore, add 30% to the dH^(z) obtained from 

Eq. (O. 

Figure la shows both GP and NPS reconstructions of the expansion history from the cosmic chronome¬ 
ter data using the best-fit values for the hyperparameters (GP) and A (NPS). Here, results correspond 
to reconstruction with the prior mean function being a constant for GP and that with the initial guess 
model being the best-fit ACDM for NPS, respectively. For completeness, we also plot H{z)/{1 + z) as & 
function of z in Figure lb. Clearly, the NPS reconstruction shows a minimum (deceleration/acceleration 
transition) at Zt ~ 0.6 whereas the GP function presents no transition in the z interval considered. This 
could be understood as the latter reconstruction might provide a more reasonable reconstruction and 
better calibration to the JLA sample if we believe that the transition from deceleration to acceleration 
has happened in the real case (also in the standard ACDM model). 


B. Calibrations for the light-curve fitting parameters 

In order to calibrate the light-curve fitting parameters and construct cosmological mo del-independent 
Hubble diagrams for the JLA SN la sample, we first transform the H(z) reconstructions above into 
distance following the approach proposed in Ref. 0. Assuming a spatially flat universe, we solve 
numerically the distance integral and calculate the corresponding uncertainty with a uniform step, Az = 
Zi+i — Zi = 0.005. For completeness, we check the difference between the results from this numerical 
treatment and the analytical calculation considering a given cosmology. At the redshift range considered, 
z < 1.2, and assuming the standard AGDM model, the difference in distance modulus is < 0.003 mag, 
which is negligible when compared to the uncertainties of current SNe la observations (~ 10“^ mag). 

Since the expansion rate measurements obtained from the derivative of redshift with respect to cosmic 
time, i.e., cosmic chronometer, and the non-parametric methods are cosmology free, distances derived 
from the reconstructed functions of Hubble parameter with respect to redshift are considered to be the 
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true ones, (luminosity distance can be obtained from the comoving distance by multiplying 

(1+z): cLl = (l+z)c?c)- And then, a and /3 are fitted using the standard minimum-x^ route to analytically 
marginalize the absolute magnitude of a fiducial SNe la, Mb, 

X^ia, 13, Mb) = A-2 X Mb X B + M^x C, (11) 

where 


A{a,l3) 

740 

= E 

a, /3; M = 0) - Mtrue(z^)J^ 

(12) 


B{a,l3) 

740 

= E 

i=l 

[/J,^^(zi, a, /3; M = 0) - fj.true(zi)} 

(13) 


B{a,P) 

740 

= E 

2=1 

1 

(14) 



Here CTjot are propagated from both the statistical uncertainties in SNe la and those in the derived 
X^ia., f3, Mb) in the Eq. (HU has a minimum at Mb = BjC [l^, and it is 

x\a,P)=A{a,P)-^^. (15) 

C{a,l3) 

Therefore, by minimizing x^(a,/3), we can achieve calibrations for a and /3 with true distances derived 
from cosmic chronometer observations. On the other hand, we also can benefit an estimation for the 
nuisance parameter, Mb, from this fitting. It should be noted that this merit enable us to break the 
degeneracy between the Hubble constant Hq and the absolute magnitude Mb when we investigate cos¬ 
mological implications of this Hubble diagram in the following analysis. 


III. RESULTS 

The light-curve fitting parameters calibrated from the reconstructions of H{z) with different prior mean 
functions for GP and initial guess models for NPS taken into consideration are summarized in Table U 
and Table im respectively. As shown in Table U it is suggested that GP reconstructions and following 
calibrations and cosmological implications are somewhat sensitive to the assumption of the prior mean 
function. However, from the Table m it is found that NPS reconstructions and following results are not 
dependent on the initial guess model. The difference between the distance modulus ^{z) derived from 
the reconstructed H{z) functions with the GP and NPS presented in Figure la, and the ACDM values 
is shown in Figure 2. In agreement with the results of Figure lb, the NPS values result in a Hubble 
diagram more consistent with the ACDM prediction than do the GP values. Here, again, it should be 
noted that, rather than fixing Hq = 70 km s“^ Mpc~^ to obtain the values of Mg and Am in the global 
fit for the ACDM model [^, these values, in our analysis, can be estimated from the cosmological- 
mo del-independent calibration by confronting light-curve fitting parameters-dependent distances of SNe 
la (Eq. [IJ with the ones derived from the reconstructed function of H{z). 

We also investigate some cosmological implications of these model-independent Hubble diagrams for 
the JLA sample. Assuming a spatially flat ACDM scenario, whose expansion rate is given by H{z) = 
77o[^m(l + -z)^ + (l —where Hq and Vim are the present value of the Hubble parameter and matter 
density parameter, respectively, we derive the constraints on the Vim — Hq plane and the results are also 
summarized in Table U and Table [Hi Once again results based on GP are slightly dependent on the prior 
mean function. In contrast, results based on NPS are hardly sensitive the initial guess model. Moreover, 
we also present constrained contours on the Vim — Hq plane with reconstructed functions in the Figure 
la taken into account. For the Gaussian processes, we obtain that the value of the Hubble constant is 
fairly compatible with the constraint from the latest PZancfc-|-WMAP9-|-BICEP2 CMB measurements 
reported in Ref. [i^. For the smoothing method, we found that the constraint on the Hubble constant 
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FIG. 2: The estimated distance modulus difference calculated with the GP, NPS and the ACDM values obtained 
from the light-curve fitting parameters. 
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is excellently consistent with the recent measurement of the local Hubble parameter obtained from the 
recession velocity of objects around us {Hq = 73.80 ±2.40 at 68.3% confidence level) [i^ (see Figure 3b). 


Prior mean functions 

Calibrations on (a, /3, M, AM) 

Constraints on Qm and Ho 

Constant 

0.137, 3.036, -19.081, -0.056 

0.29ll° °f7, 68.924t° «6 

E-D 

0.149, 3.157, -19.042, -0.029 

0.27070.838l°;g3 

ACDM 

0.158, 3.439, -18.983, -0.031 

n 9c:c:+0-029 79 '7Qr+0.521 

U.ZOtJ o.o28 5 ' ' ^^-0.522 

wCDM 

0.152, 3.307, -19.053, -0.042 

U.ZO/_o.o28 5 ''^•^-*-0-0.498 


TABLE I: Summary of the results with different prior mean functions for the GP. 


Initial guess models 

Calibrations on (a, /3, M, AM) 

Constraints on Qrn and Hq 

E-D 

0.154, 3.190, -19.022, -0.019 

0.264l°;“f7, 71.738l“:i°° 

ACDM 

0.155, 3.221, -19.014, -0.017 

0.26ll° “f7, 72.0701“:®°^ 

wCDM 

0.153, 3.172, -19.028, -0.021 

n 9«c:+0.028 y-i r9r-|-0.498 
U.ZUO_o 027’ ^ -*-•'-'^'-'-0.499 


TABLE II: Summary of the results with different initial guess models for the NPS. 


IV. CONCLUSIONS 

As is well known, implications derived from current SNe la analyses, where the light-curve fitting 
parameters are usually determined to the global fit in the frame of the standard dark energy model, are 
cosmological-model-dependent. In this paper, we have applied two non-parametric methods to reconstruct 
the Hubble expansion using 15 H(z) measurements (z < 1.2) from cosmic chronometers and transformed 


















FIG. 3: Left: Constraints on the spatially flat ACDM scenario using the model-independent Hubble dia¬ 
grams of the JLA compilation. Right: Marginalized distribution of the local Hubble parameter Hq from the 
H( 2 ;)-calibrated Hubble diagrams of JLA sample using Gaussian processes (blue solid) and the non-parametric 
smoothing method (black dashed). The dark green and dark yellow rectangles correspond to the 68.3% interval of 
Hq obtained from the latest Plancfc-I-WMAP9-I-BICEP2 CMB [3| and Cepheids observations |3| . respectively. 


the reconstructed functions into distances by carrying out a numerical integration. The choice of this 
reduced H{z) sample is based on the arguments of Refs. , which ensures that the evolution of the 

Hubble parameter reconstructed in our analysis is neither dependent on the cosmology nor on the stellar 
population model. 

By using the derived model-independent distances we have calibrated the light-curve fitting parameters 
and constructed a completely cosmological model-independent Hubble diagram for the JLA sample. 
The results suggest that the uncertainties on the light-curve fitting parameters obtained from the H{z) 
reconstructions are almost of the same order of magnitude as the ones determined in the global fit for 
the ACDM model. Therefore, we expect the constraining power of any analysis derived from these H{z)- 
calibrated Hubble diagrams to be nearly identical to the one obtained when the global fit to a given model 
is performed. It should be emphasized, however, that cosmological implications of the Hubble diagrams 
constructed from H{z) data do not suffer with cosmological model-dependence. 

However, we have shown that these diagrams and their implications depend reasonably on the method 
used to reconstruct the Hubble evolution. Furthermore, for Gaussian processes, reconstruction and 
following analysis are obviously sensitive to the prior mean function. Assuming the spatially flat ACDM 
model, we have derived constraints on the matter density parameter and Hubble constant Hq from 
the JLA sample. For the analysis that uses Gaussian processes, it is shown that constraints on model 
parameters, Dm and Hq, vary obviously when different prior mean functions are considered. In the 
reasonable and safe case with the prior mean function being constant, it is shown that the constraint on 
Hq is quite compatible with that derived from the latest Planck-|-WMAP9+BICEP2 CMB observations 
at 68.3% confidence level. When the non-parametric smoothing procedure is applied, results are hardly 
dependent on the initial guess model and consistently favor a higher value of Hq, which are in excellent 
agreement with the recent local measurement of the expansion rate from Cepheids observations. 

A final aspect worth emphasizing is that, differently from the GP reconstruction, the expansion rate 
H{z) obtained from the smoothing method behaves similarly to the one predicted by the AGDM model, 
with the deceleration/acceleration transition happening around z ~ 0.6 (Figure lb). Therefore, taking 
the standard evolution as a good description for the late-time evolution of the Universe, the NPS results 
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obtained in this analysis seem to be more consistent than those derived from the GP method. 
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